Astrophysical Effects of Scalar Dark Matter Miniclusters 
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We model the formation, evolution and astrophysical effects of dark compact Scalar Miniclusters 
("ScaMs"). These objects arise when a scalar field, with an axion-like or Higgs-like potential, un- 
dergoes a second order phase transition below the QCD scale. Such a scalar field may couple too 
weakly to the standard model to be detectable directly through particle interactions, but may still 
be detectable by gravitational effects, such as lensing and baryon accretion by large, gravitationally 
bound miniclusters. The masses of these objects are shown to be constrained by the Lya power 
spectrum to be less than ~ 10 4 Mq, but they may be as light as classical axion miniclusters, of the 
order of 1O _12 M0. We simulate the formation and nonlinear gravitational collapse of these objects 
around matter-radiation equality using an N-body code, estimate their gravitational lensing proper- 
ties, and assess the feasibility of studying them using current and future lensing experiments. Future 
MACHO-type variability surveys of many background sources can reveal either high-amplification, 
strong lensing events, or measure density profiles directly via weak-lensing variability, depending 
on ScaM parameters and survey depth. However, ScaMs, due to their low internal densities, are 
unlikely to be responsible for apparent MACHO events already detected in the Galactic halo. As 
a result, in the entire window between 1O~7M0 and 1O 2 M0 covered by the galactic scale lensing 
experiments, ScaMs may in fact compose all the dark matter. A simple estimate is made of param- 
eters that would give rise to early structure formation; in principle, early stellar collapse could be 
triggered by ScaMs as early as recombination, and significantly affect cosmic reionization. 



I. INTRODUCTION 

Despite experimental efforts to ascertain the nature of 
dark matter, its fundamental character remains a mys- 
tery. Particles motivated by supersymmetry (Weakly In- 
teracting Massive Particles (WIMPs)), extra dimensions 
(Kaluza Klein Particles), and a solution to the strong CP 
problem (the pseudoscalar axion), are among the best 
motivated from a particle physics standpoint (see [l], Q 
for reviews). Experimental efforts searching for such el- 
ementary particle dark matter have focused on utilizing 
their strong or electroweak interactions to detect them, 
either directly through their interactions with ordinary 
matter, or indirectly through their annihilations to pho- 
tons. 

On another front, gravitational lensing has already 
proven an effective way to probe the nature of dark 
matter experimentally. Photometric monitoring of many 
stars has been used to search for gravitational lensing 
by lumps of dark matter, generically known as Massive 
Compact Halo Objects (MACHOs) % H i, H . Dark mat- 
ter dominated by compact objects such as baryonic plan- 
ets, stellar remnants or black holes is now only allowed for 
m H < 1CT 7 M Q and m H > 3OM Q. However, there re- 
mains the "microlensing puzzle" , which is that MACHOs 
are observed and appear to contribute an optical depth 
toward the Large Magellenic Cloud (LMC) too large by 
a factor ~ 5 to be accounted for by simple models of 
the stellar population. Although the results are not con- 
sistent with all of the dark matter being in the form of 
MACHOs, their result [3| is consistent with an object of 
mass M ~ O.4M accounting for 20% of the halo mass. 
This signal is also at variance with the EROS experi- 
ment results @, which exclude 0.4M Q MACHOs from 
composing more than 7% of the halo. 



Searches for MACHO-like objects in the halo continue 
d, [H, and in the future, lensing searches for objects on 
cosmological scales will be increasingly sensitive to ob- 
jects in a wider range of masses (see [I, [13 for proposals 
on cosmolo gica l scale lensing using gamma-ray bursts, 
and [111, H% Il3t [Tij for gravitational lensing from distant 
quasars). In addition, a large-aperture, wide-field space- 
borne telescope is potentially capable of monitoring in- 
dividual stars for lensing effects in a program similar to 
MACHO, but in galaxy halos orders of magnitude farther 
away than the current galactic lensing experiments (see 
[lH for lensing observations towards M87 in the Virgo 
cluster using the Hubble Telescope) . In more distant ha- 
los, the Einstein radius for lensing by a given mass is con- 
siderably larger; therefore a wider range of objects can 
produce observable microlensing, and the probability of 
events increases. 

The continued increasing capabilities in lensing exper- 
iments opens the possibility of detecting lensing vari- 
ability from a wider class of astrophysical objects, be- 
yond baryonic planets, stellar remnants and black holes, 
and has the potential for studying certain types of non- 
baryonic dark matter in detail. In particular, in this 
paper we study the effects of a type of non-baryonic, el- 
ementary particle dark matter that naturally forms into 
large, self-gravitating clumps in the early universe: scalar 
dark matter miniclusters. These structures originate 
from order unity isocurvature matter density fluctuations 
created during a second order phase transition of a very 
light scalar field, sometime after the QCD phase transi- 
tion. The miniclusters form when these fluctuations sub- 
sequently collapse gravitationally at temperatures near 
matter-radiation equality. 

It has been known for some time that the QCD phase 
transition gives rise to dense axion configurations, orig- 
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inally called miniclusters [la ], or alternatively axitons 
fl7l with masses M mc ~ 10~ 12 M Q , which are de- 
tectable by pico- or femto- lensing experiments [HJ . The 
effects of large-scale modulation of axion density was also 
studied in [2(| HH, [HJ . We consider here a similar mech- 
anism for a much wider class of theories, arising from 
a second order phase transition in either a Higgs-like or 
axion-like system, and with an associated phase transi- 
tion temperature possibly much lower, perhaps even well 
below the QCD temperature. For this class of theories, 
the gravitational effects of the miniclusters arc the main 
distinguishing effects of the character of the dark matter 
and the main experimental constraint on their parame- 
ters. 

These scalar miniclusters (which we designate ScaMs 
for short) may be in a mass range interesting for lensing 
experiments. Although their densities are typically too 
low to be detectable by current generation galactic lens- 
ing experiments (and as a result, the MACHO dark mat- 
ter constraints do not apply), they may be seen through 
longer baseline galactic lensing or cosmological scale lens- 
ing. As shown below, masses as large as ~ 10 4 M Q are 
currently allowed without contradicting constraints on 
the power spectrum from the Lya forest data. Other 
constraints on MACHO-like objects from tidal effects in 
halo wide binaries (as in [23[) also do not apply to ScaMs 
on account of the ScaMs' low internal densities. A sig- 
nificant fraction (as much as half) of the dark matter 
collapses into the compact objects initially, so it is natu- 
ral to find significant microlensing rates under the right 
conditions. Unlike true MACHOs however, these objects 
are not point-like gravitational sources, but are extended 
objects with a radius which can be on the same order 
or larger than the Einstein radius, depending on the dis- 
tance to the lens. This leads to the possibility of unique 
and distinctive gravitational lensing signatures. In the 
strong lensing regime, the ScaMs produce classical caus- 
tic events with sudden appearance and disappearance of 
images, associated with sudden large-amplitude varia- 
tions in image brightness in variability surveys. In the 
weak lensing regime, sources passing behind a ScaM ex- 
perience variable small-amplitude modulation depending 
on the projected surface density of the dark matter. 

There are a wide range of theories which generate such 
objects. The second-order scalar field phase transition 
was introduced in ref. pij within the context of axion 
cosmology to show that in the presence of late phase 
transitions the axion mass and coupling constant may lie 
outside the window prescribed by the conventional as- 
trophysical and cosmological constraints; in particular a 
string scale axion with decay constant / = 10 16 GeV is 
allowed. In this model with a Higgs-like potential with a 
very small mass term, the field remains in the unbroken 
phase until the curvature of the potential is sufficiently 
large to overcome the Hubble friction, at which time the 
symmetry breaks as the field evolves to the minimum of 
the potential. Pseudo-Goldstone bosons have also been 
invoked for a variety of uses, the most famous of which 



is the axion to solve the strong CP problem. They have 
been used in relation to attempts solve the cosmological 
constant prob lem (2f|, explain the origin of large-scale 
structu re l26l |27| , and provide a warm dark matter can- 
didate [28tj . The considerable increase in our knowledge 
of cosmological parameters has ruled out or disfavored a 
number of these scenarios, however the pseudo-goldstone 
boson remains a viable dark matter candidate. 



The microlensing experiments have the capability of 
detecting these scalar objects which may never be ob- 
served directly through particle interactions due to their 
very weak couplings. The weak couplings are required 
by naturalness arguments. If there is no symmetry to 
protect their masses, radiative corrections tend to force 
scalars to be as heavy as the cut-off scale, which, for elec- 
troweak SUSY breaking, is m susy ~ 100 GeV. In order 
to maintain a scalar as light as Aq CD /M p i ~ 10~ 21 GeV, 
this scalar must be protected from this SUSY breaking; 
this is done by requiring a small enough coupling to the 
visible sector that sufficiently small radiative corrections 
are generated for the light scalar field. We require in 
particular that its coupling A to all ordinary matter sat- 
isfy A < m^/rrisusy. Such an object is truly dark, unde- 
tectable even by the most sensitive particle dark matter 
detectors; like black holes, we can detect their presence 
only through their gravitational interactions. 



In addition to lensing, these objects seed non-standard 
bottom-up hierarchical structure formation. The phase 
transition injects a large amount of fluctuation power on 
small scales so that nonlinear clusters are predicted al- 
ready at recombination; they can accrete baryons, and 
potentially trigger star formation (possibly assisting early 
re- ionization) , much earlier than the standard dark mat- 
ter with only inflationary perturbations. This also im- 
plies that measurements of the matter power spectrum 
from the Lyman-Alpha forest will limit the masses of 
these objects; we derive constraints below. 



The plan of the paper is as follows. In the next section 
we describe the physical dark matter models of interest, 
and derive general expressions for the density fluctua- 
tions. We then discuss the phenomenology of ScaMs, 
including the limits on their masses and radii from mi- 
crolensing experiments and measurements of the power 
spectrum of the Lyman-Alpha forest. We use an N-body 
code to determine their density profiles, and consider 
their gravitational lensing effects in two regimes, strong 
and weak lensing, corresponding to distant and nearby 
halos respectively. We briefly discuss the evolution of 
these objects and their impact on structure formation, in 
particular early star formation. We close by surveying 
the viability of detecting and studying this dark matter 
candidate through its lensing effects. 
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II. SCAMS AND A LATE COSMOLOGICAL 
PHASE TRANSITION 

A. Formation of isocurvature fluctuations 

Consider the generic potential of a complex scalar field 



pe 



*/f 



V{p, a) = -m 2 p p 2 + Ap 4 + /i 4 1 



/ 



(1) 



Such a potential is similar to that generated for the QCD 
axion, the difference being that QCD instanton effects 
result in /i effectively being a time dependent quantity. 
That such a potential generates 0(1) fluctuations in the 
resulting dark matter has been known for some time [l(| • 
We review the basics here. If the Peccei-Quinn sym- 
metry breaks after inflation has already occurred, it is 
expected that the initial value of the field will vary spa- 
tially, — 7r < 8 < ir. As the QCD phase transition nears 
at T ~ 1 GeV the axion gains a mass due to QCD instan- 
ton effects, generating a potential for the axion favoring 
8 — 0. Initial spatial fluctuations in the field, 8t, are then 
translated into spatial fluctuations in dark matter den- 
sity, pdm oc 8\. As the fluctuations in 8i are 0(1), the 
axion dark matter fluctuations are also expected to be 
0(1). Thus the matter power spectrum enters matter- 
radiation equality already non-linear. These fluctuations 
immediately collapse at matter-radiation equality into 
dense objects, ScaMs. For the particular case of axions, 
the mass of these objects was shown to be approximately 
1CT 12 M ~ p a (T ~ 1 GeV)d 3 H , the total dark matter 
mass within the horizon, dn, at the time of the QCD 
phase transition. 

These arguments carry over generally to any type of 
pseudo-Goldstone boson, provided, like the axion, the 
symmetry for the radial mode breaks after inflation, and 
the angular mode is effectively massless when this sym- 
metry breaking occurs. Like the axion, a phase transi- 
tion generates large density fluctuations, but we allow 
the transition to occur potentially much later, giving rise 
to more massive ScaMs. In the appendix we give an ex- 
ample of a specific supersymmetric model which gives 
rise to such late phase transitions. We call these pseudo- 
Goldstone boson modes axion-like modes on account of 
the similarity to the axion itself. These axion-like modes, 
however, need not be connected to QCD physics in any 
way. 

In similar fashion, we may also consider Higgs-like 
modes. If the radial mode is sufficiently light, a Higgs- 
like potential may also generate a late second-order phase 
transition when p breaks the U(l) symmetry. We will 
see that this second-order phase transition also results 
in 0(1) density fluctuations in dark matter condensate. 
The radial field p remains in an unbroken phase at the 
origin, pi = 0, until the curvature there exceeds the Hub- 
ble friction, m p > H , when the field rolls out to its true 
minimum at (p) = /. We will show that quantum fluc- 
tuations about pi result in 0(1) variations in the the 



roll-off time of the scalar field. These 0(1) variations in 
the roll-off time translate to 0(1) variations in the den- 
sity, due to spatial modulations in the temperature when 
the dark matter condensate forms and begins to redshift. 
In particular, we calculate 
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where bars denote mean values, T tra ns is the phase tran- 
sition temperature in a particular Hubble patch, and 
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It will be useful to rewrite this in terms of a time delay, 
5t from the mean time of the phase transition, ttrans, 
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where we have included the explicit spatial dependence 
of Sp(x) on scales exceeding the Hubble size at the time 
of the phase transition. We can see that in the limit that 
8t < t, the result reduces to 
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which is identical to the inflationary result for an upside- 
down harmonic oscillator [29j , up to the multiplying fac- 
tor. 

As a result, in order to calculate the density fluctu- 
ations, we need only calculate the time delay. This is 
approximated by [291 ] 



St(x) ~ lim 
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where we have expanded in fluctuations, 6<fi(x, t), around 
the background field, cp^, 0(x, t) = (o) (i) + 6<f>{x, t). 

We determine the classical evolution of the background 
field (j)^ from its equation of motion, 



4> {0) (t) + 3ff0 (o) (t) + m 2 (o) (t) = 0. 



(7) 



Since we are interested in the behavior at asymptotically 
late times, we ignore the Hubble friction terms, and the 
solution is 



(8) 



Likewise, we can determine <5</>(x, t) from an equation of 
motion 

<50(x, t) + 3H6(j)(x, t) + m 2 (50(x, t) - V 2 (50(x, t) = 0. (9) 

For large times, we may neglect the spatial gradient and 
Hubble friction terms, and the solution is 



50(x,i) = <50(x,i = 0)e" 



(10) 
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Thus St(x) is set entirely by initial conditions: 

Kx,t = o) 1 



6t(x) 



c/>( )(t = 0) m 



(11) 



We will assume that the background field <fi^ is ini- 
tially at the origin. However, we cannot choose (j)^ (t — 
0) = because of quantum fluctuations. The size of these 
fluctuations are calculated from the Fourier transform of 
the two point function, following (29|: 



A0(k) EE 



d 3 a 



_(2tt) 3 J (2tt) 3 

d 3 k' f d 3 x 



e lk - x (^(x,i = 0)50(0)) 

1/2 



1/2 



(2tt) 3 J (2tt) 



.gik-Xgik'x 



(32tt 3 ) 1 /2 



(12) 



For the mode of interest, k — m, these fluctuations set 
both (t = 0) and the root mean square (rms) value 
of 5</>(k), the Fourier transform of 6<p(x,t = 0). Putting 
it all together, we have the rms density fluctuations, 
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showing that the density fluctuations are of order unity. 



B. ScaM collapse 

These fluctuations collapse gravitationally into ScaMs 
around the time of matter-radiation equality. For col- 
lapse of a uniform sphere the final core density of a viri- 
alizcd ScaM is 111 



p = 140<5 3 (<5 + l)p eq , 



(14) 



where S = Sp/p. Their masses are set by the dark matter 
mass inside the horizon at the time of the phase transi- 
tion, 



MscaM = -^TrdffPDAliTtrans), 



(15) 



where dn = H(T trans )~ x . This corresponds to a ScaM 
mass, assuming these scalars compose a fraction r of the 
dark matter, 
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C. Limits on ScaM Mass from Lyman-a 

As a result of these non-linear density fluctuations, 
the phase transition adds a large amount of fluctua- 
tion power into the spectrum on small scales. These 



scales, r s ~ T trans /(T H (T trans )) (~ 10 pc (comoving) 
for T trans ~ A-qcd) are well below the reach of current 
measurements of the power spectrum, since the small- 
est scale measurements, derived from observations of the 
Lyman-q absorption of the spectra of distant quasars 
[30L l3lL l3^ |. reach down only to scales ~ 0.1/i _1 Mpc, 
the scale where protogalactic gas is collapsing into mildly 
nonlinear filaments. The phase transition, however, gen- 
erates a Poisson white noise power spectrum on scales 
larger than r s which adds to the inflationary power on 
the Lyman-a scale, 
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where ns C aM is the number density of ScaMs and the 
subscript p is for primordial. The added power today 
is then the product of the primordial white noise power 
spectrum with the transfer function for isocurvature fluc- 
tuations, T iso , 



where 
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Now neither P p nor Ti SO is wave number k dependent, 
whereas Pcdm decreases with k. We plot this power 
spectrum in fig. 1. We have introduced in the white 
noise spectrum a smoothing scale r s ~ du {T trans ), on 
which the Kibble mechanism smooths field fluctuations, 
to remove power on the smallest scales, 
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We can see from fig. 1 that for sufficiently large k, 
PscaM will exceed Pcdm', we must ensure that this oc- 
curs on smaller scales than are reachable with Lyman- 
alpha measurements, k > kj, so that we require 

h 1 + Zeq fM^L K P Lya ( kj ), (21) 

4 PDAI 

where kj ~ lOh Mpc -1 . This yields the constraint 

MscaM < 4 x 10 3 M Q , (22) 

corresponding to a constraint on the temperature of the 
phase transition, 



2 x lO -0 GeV < T t , 



(23) 



assuming the scalars compose all the dark matter. A 
similar constraint from the Lyman-alpha power spectrum 
was derived in [33] using numerical simulations in the 
context of primordial black holes. 

Phase transitions between the QCD scale and this 
Lyman-a limit create abundant ScaMs in the mass range 

1O- 12 M < M ScaM < 4 x 1O 3 M , (24) 

much of which is accessible to current and future mi- 
crolensing experiments. 
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FIG. 1: CDM power spectrum derived from BBKS HH with 
white noise power spectrum smoothed on scales r s ~ 7 x 
10~ 2 Mpc, given by eqn. 1201 The amplitude of the white noise 
spectrum corresponds to the power spectrum for a ScaM of 
mass M~4x 10 3 M Q . 



D. ScaMs as microlenses 

For ScaMs to act as observable strong microlenses, 
three conditions must be satisfied: 

1. ScaM masses must lie in the mass range reachable 
by experiments. For classic microlens searches by 
stellar monitoring in the local halo, this range is 
currently 10~ 7 Af o < M ScaM < 1OM . However, 
the accessible range will widen in the future as mi- 
crolensing experiments access fainter and more dis- 
tant monitored background sources in other galax- 
ies. 

2. A rough criterion for strong lensing, leading po- 
tentially to large-amplitude variations in source 
brightness, is that the radius of the ScaM be smaller 
than the Einstein ring radius. For microlensing of 
objects over cosmological distances, the Einstein 
radius is 
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For lensing toward a source in the local galactic 
neighborhood (e.g. toward the LMC or M31), the 
Einstein Radius is 
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where D is the distance to the lens, and it is as- 
sumed that D <C D s , the distance to the source. 

Using the spherical model, eqn. [141 and assuming 
that the ScaMs are approximately constant density, 
we calculate 
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3. Their cosmological abundance must be consistent 
with the limits from the current lensing experi- 
ments. Since these objects would generally be too 
fluffy to create strong lensing in the nearby halos 
observed by the current generation of galactic mi- 
crolensing experiments, consistency with the limits 
of these experiments is generally not problematic. 

Although these objects would generally not be dense 
enough to be observed by Galactic microlensing exper- 
iments, they may be detectable as microlenses for more 
distant sources and halos. If they do produce strong 
lensing events, they do not obey the classic Paczynski 
[HI point-mass light curve, but instead are dominated 
by more generic caustic-crossing events. More generally, 
in nearby halos they may not even act as strong lenses, 
but may have a resolved density structure that appears as 
small-amplitude variations in the light curve of a lensed 
source. 

To improve on the spherical collapse model, eqn. 1141 
and in particular to determine properties of these ob- 
jects observable by lensing experiments, we simulate the 
collapse of ScaMs using an N-body code. The resulting 
objects are more realistic than the spherical model and 
allow determination of some representative density pro- 
files. 



III. SIMULATING SCAMS 

We simulate the formation of ScaMs in the radiation 
dominated era using the N-body code described in [36j. 

A. Initial conditions 

The initial density profile may be determined utilizing 
one of two methods: either by solving the classical equa- 
tions of motion for a field <fi directly, or simply using the 
power spectrum of eqn. 1171 The equation of motion for 
a scalar field is given simply by 



3H(f) 
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where the Laplacian is taken with res pec t to comoving 
coordinates x. We can rewrite this (see [18| for details) in 
terms of conformal time, 77 = R/Ri, where R\ is defined 
by H(Ri) = ?7i0, and comoving Laplacian V 2 , taken with 
respect to coordinates x = H(Ri)Rix, 
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We assume the system is subject to white noise initial 
conditions, 

sinful) . 
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where the £'s are random phases. 

We solve the equation of motion on a lattice 100 
sites per side for the axion-like potential, V(4>) = 1 — 
cos(0//) = 1 — cos 6. This potential would result in the 
formation of domain walls if, for example, A — irf (in the 
language of axions, this corresponds to multiple vacua, 
N > 1). The walls are transient objects that quickly dis- 
sipate by particle radiation, but introduce singularities 
that make numerical integration difficult. Since we are 
interested in the formation and evolution of the ScaMs 
only, we avoid this problem by choosing A so that 9 varies 
between —n and 7r, and the root-mean square (rms) value 
of 9i — <pi/ f is the rms average value of the misalingment 
angle, 9 rms = ir/\/3. In this case, the domain walls never 
form in our box. 

Solving the equation of motion numerically automat- 
ically simulates the effects of the Kibble mechanism, 
smoothing field fluctuations on scales smaller than the 
horizon size at the time of the phase transition. We plot 
in fig. [2] a two dimensional slice of the initial white noise 
density fluctuations for the axion-like potential. 

We evolve these fluctuations to 77 = 10, at which point 
the fluctuations are expected to remain mostly spatially 
frozen (modulo logarithmic growth of fluctuations) un- 
til gravitational collapse begins right around (or even 
somewhat before for the most dense ScaMs) the epoch 
of matter-radiation equality. The density fluctuations in 
a box with sides whose length are four times the hori- 
zon size at the time of the phase transition are shown in 
fig. O The density fluctuations have been normalized to 
the average density in the box, so that p(x)/p is shown. 

We choose an alternate method to determine initial 
density fluctuations for the Higgs-like potential. Ac- 
cording to eqn. [T31 the Higgs-like potential gives only 
quasi-nonlinear density perturbations, which allows us to 
model the formation of these clumps realizing the initial 
density fluctuations using the standard N-body particle 
method: particles are initially placed on a lattice to min- 
imize shot noise, then displaced from those positions by 
adding perturbations, mode by mode in the Zcl'dovich 
approximation, with amplitudes and phases selected ac- 
cording to the distribution derived from the power spec- 
trum. We adopt the power spectrum of eqn. [^D] 



P(k) = Ae 



-(kr s /2ir) 2 /2 



Mpc d 
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which creates white noise filtered on a scale r s . This 
smoothing scale corresponds roughly to the horizon size 
at the time of the phase transition. The expected density 
fluctuations on this scale are thus 




FIG. 2: Two-dimensional slice in the x-y plane of initial white 
noise field energy density distribution, x-y coordinates are 
in 77, where 77 = 1 corresponds to a length d,H(Tt r ansl) one 
horizon size at the time of the phase transition. The z-axis 
is the initial white noise over-density p(x)/p, where p is the 
mean density in the box. 




FIG. 3: Density distribution after the phase transition, at 
77 = 10 for the axion-like potential. Axes are same as in 
fig- [21 Note the highly nonlinear nature of the initial den- 
sity perturbations. This distribution remains fixed until near 
matter-radiation equality when it evolves gravitationally into 
collapsed ScaMs; this density spectrum is the input for the 
N-body simulation. 
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where k s = 2Tr/r s and we choose r s so that the result- 
ing ScaMs haves masses around 1M Q (calculated from 
MscaM ~ Pdm^s, and corresponding to a phase transi- 
tion temperature, T tra ns ~ 10~ 4 GeV), and A such that 



dk 



(5p/p) rms — 0.5. Note that while we use specific phys- 
ical scales for the purpose of the simulation, the result 
is expected to be completely scale invariant, and should 
apply with suitable rescaling to ScaMs of all masses. 

We plot the corresponding density fluctuations for the 
Higgs-like potential in fig. 0J We then use the N-body 
code to evolve the objects into the collapse epoch near 
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FIG. 4: Density distribution generated by the smoothed white 
noise power spectrum. The size of the density fluctuations 
were chosen consistent with the analytic result for the Higgs- 
like potential of Sec. II. Like fig. [3] this distribution is input 
for the N-body simulation. 



matter-radiation equality. 



B. N-body simulation of ScaM density profiles 

Using the density profile generated by solving the clas- 
sical equations of motion for the axion-like potential, a 
grid of 100 3 particles is laid with masses weighted accord- 
ing to the locally computed initial density. The evolution 
of the particles is started at a redshift 1 x 10 5 ; objects 
are fully collapsed by matter radiation equality, around 
a redshift of 3000. As the density fluctuations are ini- 
tially much smaller for the Higgs-like distribution shown 
in fig. U the evolution can be started much later; we 
choose a redshift of 1 x 10 4 , and objects are fully col- 
lapsed by a redshift of 1000. 

We show in figs. [51 El a slice of the final particle distribu- 
tions of the collapsed objects, at redshifts around matter 
radiation equality. The initial over-densities shown in 
figs. O [H evolve into the fully collapsed objects shown in 
figs. 03 [H The pictures of course resemble images of the 
low-redshift cosmic web of dark matter, but in this case 
represent very small-scale, nonlinear objects at a redshift 
around 10 3 . The "normal" inflationary perturbations in 
dark matter at this time still have a very small amplitude 
on all scales, of the order of one percent. 

We plot in figs. [71 [5] the density profiles of a sampling 
of ScaMs, where we choose to normalize our densities 
and radii against those predicted by the spherical model, 
eqn. Qj] for 6=1. We can see that the naive spherical 
model is approximately accurate in predicting the maxi- 
mum density of the ScaMs. The axion-like potential has 
typical overdensities 6 ~ 5— 10, as shown in fig. [31 so that 




FIG. 5: Snapshot of structure formation at z = 3000 for 
axion-like ScaMs; the plot is colored according to the loga- 
rithm of the density. The scale here is 1.2/i _1 Mpc, in co- 
moving coordinates, although the result is expected to be in- 
variant for any box size, given that we rescale the smoothing 
scale accordingly. Note that structure formation has already 
commenced before matter radiation equality. 




FIG. 6: Same as fig. [5] but for the Higgs-like potential at a 
redshift z = 1000. This simulation was evolved to a lower 
redshift, as the initial over-density of the Higgs-like system is 
lower, and hence gravitational collapse occurs later. Although 
difficult to see from this rendering, the densities of the ScaMs 
are also lower. 



the spherical model predicts p/p sp h(S = 1) ~ 100 — 5000. 
The Higgs-like potential, on the other hand, has typi- 
cal overdensities S ~ 0.5 — 1 (fig. [3}, and so we expect 
p/psph(S = 1) ~ 0.1 — 1, again consistent with the cen- 
tral densities given by the simulation. The ScaMs with 
higher central densities correspond to ScaMs which had 
higher initial overdensity 5. 
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-2 I — 1 1 — i — — . JE^, 

-1.75 -1.5 -1.25 -1 -0.75 -0.5 -0.25 
log(r/R Bph ) 



FIG. 7: Axion-like ScaM density profiles for five gravitation- 
ally collapsed ScaMs. Vertical axis is \og(p(r)/p sp h), where 
Ps P h(S — 1) — 280/5 e? ; that is, we normalized the density 
profile against the uniform density prediction of the spherical 
model, eqn. 1141 with S — 1. Horizontal axis is \og(r/R sp h), 
where R 3p h is the radius computed from p sp h and the total 
mass contained within the horizon at the phase transition, 




-1 -0.8 -0.6 -0.4 -0.2 

logfr/H.pi) 



FIG. 8: Same as fig. [7J but for the Higgs-like potential. 

C. Strong lensing profiles 

What is the implication for lensing experiments? As 
explained in section III, in order for an object to act as 
a strong gravitational lens, the enclosed mass, M enc i, at 
any given radius must exceed a minimum, 

M encl (r) >lM & ^-^j, (33) 

where s(D) is dependent on the base length for lensing, 
D (s - 3 x 10 15 cm for D = 5 Mpc, and s ~ 3 x 10 14 cm 
for cosmological scale lensing) . Equivalcntly, at any given 
distance from the center of the object r, the radius must 
not exceed the Einstein radius, 

The Einstein radius can be computed directly from den- 
sity profiles shown in figs. [5) 



We plot in figs. [H [10] Re versus r, for s — 3 x 
10 15 cm , 3 x 10 14 cm. As in figs[7J [SJ we have normal- 
ized the radii against R sp h, the prediction of the spherical 
model. The objects lens if Re > r, that is, if the Einstein 
profiles lie above the straight line, Re — t, shown in the 
figure. Although the simulation was run for a particular 
mass within the horizon at the time of the phase transi- 
tion, M ~ 1M©, it is simple to apply the results to many 
different ScaM masses simply by rescaling the total mass 
in the box. In this scaling Re increases like M 1 / 2 , but 
the box size decreases like M 1 / 3 so the net result is that 
the vertical axis in the figure is scaled up by M 1 / 6 , mean- 
ing that more massive ScaMs lens more easily. Note that 
the mass labeled is the mass enclosed within the horizon 
at the time of the phase transition, not the lensed mass 
of the ScaM, which may be significantly lower. 

We can see that both Higgs-like and axion-like ScaMs 
are of interest for lensing experiments on a cosmological 
scale, from objects as light as 1O~ 6 M0 to as heavy as 
the cosmological bound of 4 x 10 3 M Q . While neither of 
these objects is of interest for the current galactic lens- 
ing experiments with baseline D ~ 50 kpc, the axion-like 
configuration in particular will be of interest for lensing 
in the Mpc range and greater, as might be accessible in 
future experiments. We also note here that we have not 
chosen the most extreme set of initial conditions to sim- 
ulate for the axion-like potential. The axion itself has 
a time dependent mass around the time of the phase 
transition which results in larger density contrast. Nor- 
malizing against the size of the box, we find that axions 
in particular generate overdensities S ~ 20 — 30. (Kolb 
and Tkachev [18| found larger overdensities 6 ~ 30 — 100 
as they normalized against the density for 9 rms — tt /\/3 
and not the actual average density in the box.) In this 
case, Re for the axion-like potential is increased by as 
much as a factor of 10, which improves the observational 
prospects. 

As long as the radius of the lens is much bigger than 
the size of the background source, those situations where 
Re > r will occasionally lead to strong-lensing events 
with high amplification. The frequency of this happen- 
ing is roughly given by the fraction of the halo mass in 
ScaMs above this threshold, times the mean lensing op- 
tical depth of the halo, times A~ 2 where A is the amplifi- 
cation. The mean optical depth is a very small for nearby 
halos (ss 10 -6 for the LMC experiments) but is typically 
« 0.1 for halos at cosmological distances. Such consid- 
erations affect the number of sources and the cadence of 
observations needed in a variability survey. 



D. Weak lensing: mapping ScaMs using variability 

Even if ScaMs are not dense enough to cause strong- 
lensing events, they can in general still produce weak 
lens amplification. In some situations the systems are 
small enough that this amplification is time variable, and 
the projected density profiles of ScaMs can be measured 
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FIG. 9: Axion-like ScaM Einstein radius for the enclosed 
mass, derived from the computed N-body density profile, ver- 
sus ScaM radius, again normalized against the ScaM radius 
Rsph predicted by the spherical model. The diagonal line di- 
vides where Re > r, when lensing is possible. The upper set 
of curves in each plot is the Einstein radius for cosmological 
scale lensing, the lower set of curves the Einstein radius for 
lensing at a distance D = 5 Mpc, which may be relevant for 
MACHO-type experiments in nearby galaxy halos. This is 
shown for three different ScaM masses, marked at the lower 
right in each plot. 



directly by monitoring sources. Survey parameters, such 
as the distance of the halo under study and the number of 
sources to be monitored, can be optimized depending on 
the predicted ScaM parameters; in principle, both weak 
and strong lensing can be studied at various distances. 

Consider a standard gravitational lens mapping a 
source at position 9s in the source plane to a position 
9 j in the image plane. In the image plane the mapping 
is characterized by the convergence n and complex shear 
7. The magnification fi is the inverse Jacobian of the 



FIG. 10: Same as fig. [9] but for the Higgs-like potential, 
mapping [37| 

M" 1 = (1 - k) 2 - | 7 | 2 . (35) 

The convergence is determined by the surface density £ 
along the line of sight, 

K = E/E C) (36) 

where the critical surface density depends on the (angular 
diameter) distances to the source (Ds), the lens (Dl), 
and between them (-Dls) as: 

S c = D ls c 2 /AtiGD l D s « 3.5kg m- 2 {D LS Gpc/D L D S ). 

(37) 

If £ > Sc, generally light rays focus to a point some- 
where, the mapping involves multiple images and gener- 
ally strong shear, and amplification is both nonlocal and 
nonlinear. High-amplitude variability is in the strong 
lensing regime and dominated by caustics near fold catas- 
trophes in the mapping. By contrast, in the case of weak 
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lensing, S < £c, the amplification of a source is more di- 
rectly related to the surface density in that direction. For 
small surface densities, and where the shear is negligible, 
the magnification is just 



d8 s 



86, 



1 + 2k. 



(38) 



Thus as a source traverses behind a ScaM and X varies, 
its brightness changes by a fraction 2S/Sc- We plot in 
figs. QTJ Q2] S/Ec for the axion-like and Higgs-like po- 
tentials. We see that the variability will be difficult to 
detect for the current lensing experiments towards the 
LMC and M31, as the effect is only a few percent even 
for the most massive ScaMs. The weak lensing variabil- 
ity, however, will be of interest for longer baseline ex- 
periments, D = 5 Mpc, and also for cosmological scale 
lensing. For an axion-like ScaM of mass 1M Q , the time 
scale for variability is on the order of a year (assuming 
a ScaM velocity across the source of 300 km/s), and de- 

1/2 

creases with ScaM mass as M ScaM . 

Note that the weak lens effects cover more area, and 
have a larger probability of affecting a background source 
than the strong lens effects. For a given mass of lensing 
material, the lensing probability or optical depth scales 
in the weak regime like r oc oc (/i — 1) , but the 
distribution of mass in the outer parts of ScaMs is even 
more favorable to the weak lensing program. In the ex- 
amples shown in figures 11 and 12, we see that variability 
at the ten percent level (that is, S/Sc ~ 0.05) occurs at 
a radius which is typically more than ten times the ra- 
dius of strong lensing, so variation of this magnitude is 
over 100 times more frequent than strong lensing events. 
Photometric accuracy and stability are the main practi- 
cal limits in mounting variability surveys around this ef- 
fect, but long term variations at the few percent level are 
certainly within the range of proposed instruments (such 
as JDEM) designed to monitor distant supernovae. 



IV. EARLY BARYON COLLAPSE 

With WIMP-type dark matter such as neutralinos, in- 
flationary fluctuations lead to collapse of the first dark 
matter halos of roughly an earth mass scale, determined 
by the damping scale of the WIMP streaming motions 
[38l |. The collapse begins on after a redshift of 100 and 
on small scales has little effect on the baryons, since the 
gravitational potentials are so shallow. 

With ScaMs on the other hand, the addition of isocur- 
vature fluctuations creates deeper dark matter potentials 
at earlier times. We have already seen that the accretion 
of absorbing gas, leading to Lyev absorption, is indeed the 
main astrophysical constraint at present on ScaM mass. 
In this section, we estimate the nonlinear effects of ScaMs 
on early structure, using an analytic Press-Schechter ap- 
proach to bottom-up hierarchical clustering. 




-1.2 



-0.8 -0.6 -0.4 

log (r/Raph > 



-0.2 



FIG. 11: Log mean surface density E relative to the critical 
density E c as a function of log radius, averaged in annuli, 
normalized again to the prediction of the spherical model, 
Rsph- E/E c ~ 1 corresponds roughly to the onset of strong 
lensing, and 2E/E C corresponds to the fractional weak- lens 
amplification. Higgs-like potential. 




FIG. 12: Same as fig. 1111 but for the axion-like potential. 
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Once the ScaMs form at z eq w 3300, bottom-up hi- 
erarchical structure formation begins. In the spherical 
model, fluctuations on a scale M collapse into structures 
once the fluctuation has grown to a size 

8{M) ~ 1.69. (39) 

As the fluctuations are isocurvature, they grow as 



onto a structure of mass 
1M 



M 



5{MY \1 + Za 



1 + z 



< <l 



(47) 



6{M) = S{Ms^)yj^\^, (40) 

where 5(Msco,m) — 1- Therefore a structure of mass M 
will form at a redshift Zf 



1 + ^"^M)V^-2 (1 + ^ ) - (41) 

By making use of the spherical approximation, pM = 
U0S{M)l(5(M)i + l)p eq , where 



S(M)i = S(M ScaM ] 



M s 



M 



(42) 



we also determine the virial radius of a structure of mass 
M, 



R(M) 



M 



1/3 



M 



(43) 



^6 X l0 2 S(M ScaM )p eq J V MscaM 

and hence the virial velocity of these structures is 
jGMjR ~ G'/ 2 (6x 10 2 S(M ScaM )p eq ) 1/6 (MM 3 ScaM )^ 12 



5 x 10" 



M 
1Mb 



1/12 



/Ms 



cnM 



V 1M G 



1/4 



(44) 



If the baryons have a thermal velocity greater than the 
virial velocity, they are prevented by pressure gradients 
from collapsing into the ScaM. To trigger accretion and 
collapse we require the ScaM to have a virial temperature 
greater than that of the microwave background, since the 
baryons can certainly not radiate binding energy if they 
are cooler than that. The velocity of the baryons at the 
CMB temperature today would be v ~ J 3To /m p ~ 9 x 
10~ 7 . The baryons begin accretion onto the ScaM seeded 
structures once the velocity of the baryons drops below 
the structures' virial velocity, which occurs at redshift, 



2 x 



1 n — 6 ( MscaM 

iU 1 1M 



1/3 



(1+z 



2S(M) {1+z, 



1/6 



(45) 



3T (l+z acc ) 



Solving this, we find that baryon accretion onto the struc- 
tures formed hierarchically from a ScaM seed of mass 
MscaM, occurs at a redshift 



l + z a 



30 



M s 



caM 



1/2 



(46) 



We do not pursue further here the physics of the ac- 
creted baryons, which depends on complex details of non- 
linear collapse and cooling [3^, S3, IHJ • However it is 
clear that the early star formation in the ScaM models 
can start much earlier than in standard CDM cosmol- 
ogy where significant collapses occur at a redshift less 
than about 30. With massive ScaMs, collapse can oc- 
cur as early as recombination. Since cooling is relatively 
efficient at such early times (for example, through Comp- 
ton cooling on the microwave background) , these models 
also likely produce very early stars. These effects may 
eventually be observable either directly via deep infrared 
imaging, or indirectly via reionization effects. In any case 
it is clear that there would be significant modifications to 
early star formation with standard CDM perturbations. 



V. CONCLUSION 

We have not attempted to trace the evolution of the 
ScaMs in detail to the present, where they would be in- 
corporated hierarchically into the standard galaxy-size 
dark matter halos. Although they would be subject to 
some disruption in the course of hierarchical assembly, 
there is good reason to think that they would mostly 
survive intact to the present. Certainly they fare bet- 
ter than the much more diffuse earth-mass halos of neu- 
tralino CDM [38| . In that case, the very flat distribu- 
tion of density leads to a large range of masses collaps- 
ing almost at the same time on top of each other. Early 
ScaMs cluster instead with a steep white-noise spectrum, 
where the mass scale grows as a power law. In that sit- 
uation there is more room for survival and less chance 
for disruption; the process of early ScaM clustering has a 
similar spectrum to, and therefore resembles, the larger 
galaxy-scale hierarchy today, where simulations have es- 
tablished that much of the satellite substructure survives. 
The dumpiness of dark matter in the halo may also have 
observable consequences through tidal forces, such as dis- 
ruption of globular clusters; these effects have been stud- 
ied in the context of black hole dark matter and other 
highly compact objects [il. lij. liEj . 

We have shown that in scalar theories with a cosmolog- 
ical phase transition below the QCD scale, ScaMs form 
with a mass in an interesting range for microlensing ex- 
periments, 10" 12 M Q < MscaM < 10 4 M Q . The radii and 
density profiles of the ScaMs vary considerably, depend- 
ing on the initial size of the density fluctuations. For an 
overdensity 6 ~ 50 and mass M ~ IMq, these objects 
would be visible for current galactic scale lensing experi- 
ments. Objects with (more plausible) lower overdensitics 
S ~ 5 — 10, while not detectable with the current gener- 
ation of experiments, would be visible for future lensing 
experiments with longer baselines, D ~ 5 Mpc or greater, 
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requiring a space based platform (or wide field adaptive 
optics) to conduct a lens-induced- variability survey. Cos- 
mological scale lensing is the most powerful tool of all 
for detecting these objects, but on that scale, detailed 
monitoring studies are still in the distant future. Some 
information might be obtained sooner from quasar mi- 
crolensing, but ScaMs tend to be not much larger than 
quasar emitting regions, so this technique has limited 
application and dynamic range. Smaller sources such as 
individual stars are, of course, much fainter. 

The scalar field composing these objects is so weakly 
coupled to the standard model that they would never be 
detected directly, but would only make their effect known 
to us gravitationally. Thus the rich substructure of the 
halo dark matter does not have other effects, such as 
gamma ray annihilation signatures or clustering in direct 
detection experiments, that occur for other dark matter 
candidates. The main exception to note is the possibility 
of ScaMs of classical invisible axions, which could show 
up in direct detection as clumping in time as well as 
energy. 

ScaM particles when they condense are moving with 
modestly relativistic velocities in spatially coherent 
streaming flows on the scale of the nonlinear lumps at 
that time. These redshift by the time of the nonlinear col- 
lapse so that theydo not prevent collapse, but it is worth 
mentioning that they do not have the dynamically-cold 
classical distribution function characteristic of homoge- 
neous axion condensation. The broader coarse grained 
distribution function could in principle have some dy- 
namical effects, which are too subtle to model in our 
simulations. 

The scalar field collapses into ScaMs at matter- 
radiation equality, seeding early bottom up hierarchical 
structure formation as successively larger mass scales be- 
come virialized. This in turn results in early star forma- 
tion, with baryon accretion onto the Scams starting as 
early at z ~ 1400 for a ScaM mass near the cosmological 
limit at Ms C aM ~4x 10 3 Af Q . This could have substan- 
tial observational effects on energy input into diffuse gas, 
affecting the epoch of rcionization, and possibly also di- 
rect detection of early stellar systems. The reionization 
epoch is already being probed by quasar absorption to 
z > 6 and by CMB polarization to z > 10; even larger 
redshifts may become accessible in the future to similar 
techniques, as well as direct 21cm mapping. These more 
indirect effects are of course more complicated to model 
than the more direct influence of the lensing. 

The nature of the dark matter remains unknown to 
us. Scalar fields play an important role in cosmology, 
being front and center in viable models of inflation and 
many theories of the dark energy. We have shown here 
that, should a scalar field also contribute to solving the 
dark matter mystery, it is a good candidate for large 
isocurvature fluctuations which cause them to collapse 
into dense ScaMs; in fact, though the scalar may live in a 
hidden sector extremely weakly coupled to us, so that we 
would detect no direct particle interactions, gravitational 



lensing may provide us a probe into this sector. 
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APPENDIX: LATE PHASE TRANSITIONS IN 
AXION COSMOLOGY 

We demonstrate a technically natural model which 
provides a phase transition well below the QCD scale, 
seeding density fluctuations in both Higgs-like and axion- 
like fields. This model was motivated in ref. [24[ for the 
purpose of demonstrating that in theories where the ax- 
ion's decay constant evolves after the QCD phase tran- 
sition, the cosmological upper bound on the axion decay 
constant, f a < 10 12 GeV, does not apply. The model as 
written here applies specifically to the QCD axion, but 
the model can be generalized to any set of scalar fields, 
not necessarily connected with the solution of the strong 
CP problem. 

The phase transition in this model is driven by a com- 
plex field <fr, whose radial mode has a Higgs-like potential, 



v(4>) = / 



M 2 



(A.l) 



where we require /j, 2 /M -C H(Tqcd) in order that the 
vacuum expectation value (vev) of the field remain frozen 
at its initial position, away from the minimum of the 
potential at (0) = M, until T < T QCD . The field (j) 
is coupled to the Peccei-Quinn (PQ) sector through a 
potential 



V(4>,X 1 ,X 2 ) = X 2 \hXiX 2 



; 2 1 2 



(A.2) 



where X% and X 2 carry opposite PQ charges. X\ is the 
standard PQ field coupled to gluon field strength through 
the term 



X\GG, 



(A.3) 



where the usual PQ potential fixes the vev |Xi| = fi 
when the PQ symmetry breaks at temperatures T S> 
Tqcd- At the QCD phase transition, the axion gains a 
mass m a ~ m^f^/ f\, and a dark matter condensation of 
axions forms. 

A phase transition occurs when ~ H(T), so that 
<f> rolls out to the minimum of its potential at \<f>\ = M; 
\X 2 \ then follows the flat direction in eqn. IA.2I to a vev 
f 2 = M 2 /hfi, where the hierarchy f 2 > M > fi is 
assumed. It can then be shown (see 2J] for details) that 
the energy density in the axion dark matter is transferred 
into the heavier pseudo-scalar ni = fiArg(Xi), having 
mass m ni AM 2 //i generated at the phase transition. 
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As a result of the phase transition, there are two dark 
matter candidates in this model, one Higgs-like and one 
axion-like. At the phase transition, the Higgs-like mode 
\<f)\ releases energy fj, 4 which subsequently redshifts like 
dark matter, and may be cosmologically abundant. It is 
so weakly coupled to ordinary matter that it may only 
be detected gravitationally. It has 0(1) density perturba- 
tions resulting from the phase transition at a temperature 
T <C Tqcd, which collapse into Higgs-like ScaMs with 
mass much larger than an axion minicluster, Mg ca M 3> 
1CP 12 Mq. The axion-like pseudoscalar tti may also be 
cosmologically populated, and collapses into axion-like 
ScaMs. The mass of these ScaMs is much larger than 
QCD axion miniclusters as the pseudo-Goldstone boson 
in this model has an additional mass generated at the 
phase transition m- Kl /m a ^~ 12 MQ 3> 10 -12 M Q . 

There is a mechanism in this model which allows for 
equal cosmological abundances of Higgs-like and axion- 
like dark matter, <j) and 7Ti, without any fine-tuning. At 
the QCD phase transition, tti dark matter is produced 
with energy density 

P-K t (T) = m^n^, (A. 4) 

where n ni is the number density, which dilutes as T 3 , and 
m 7ri now receives a contribution from its mixing with </>: 

2 mlfl A 2 (|0|) 4 

"V W 72— + 72 > ( A - 5 ) 

Jl Jl 



This creates an effective potential for 0, 



v eff(4>) = m OT1 (?!>)n 7ri . (A. 6) 



This effective potential will delay the temperature of the 
phase transition, when <j> rolls to (\4>\) = M, until 7Ti 
has diluted enough that V e ff(M) < fi A . At that point 
the energy densities in the 4> and ~k\ fields are equal, and 
they redshift concomitantly so that their energy densities 
remain equal thereafter. This gives rise to the possibil- 
ity of the cosmological presence of both Higgs-like and 
axion-like miniclusters. As the Higgs-like and axion-like 
miniclusters have different typical densities (the latter of- 
ten being 10 4 times more dense, as explained in Sec. Ill), 
this implies a potentially varied cosmic population reach- 
able by lensing experiments on different scales. 

We also wish to emphasize that the physical features 
discussed here of late scalar field phase transitions are 
quite generic and independent of the presence of axions; 
however, the axion model provides a motivation to con- 
sider such phase transitions around the QCD time, as 
well as illustrates a model with all the physical features 
discussed in this paper. 
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